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While  sampling  at  a  Nyquist  frequency  equal  to  the  highest  frequency  present  in  the  data 
(critical  sampling)  is  sufficient  to  prevent  aliasing  in  both  the  data  and  the  autocorrelation  of  a 
bandlimited  energy  signal,  the  sampling  requirements  for  the  avoidance  of  aliasing  in 
higher-order  correlations  and  spectra  are  not  the  same.  Also,  there  is  a  difference  in  aliasing 
effects  depending  on  whether  one  samples  the  original  continuous-time  signal  and  calculates  the 
autocorrelation  or  one  samples  the  continuous-time  autocorrelation.  This  distinction  between 
sampling  procedures  must  be  made  for  correlations  of  higher  order,  as  well,  for  which  not  only 
the  type  of  aliasing  but  also  the  sampling  requirements  to  prevent  aliasing  differ.  In  particular, 
if  one  samples  the  continuous-time  autobicorrelation  or  autotricorrelation,  critical  sampling  is 
sufficient  to  prevent  aliasing.  In  practice,  however,  it  is  not  usually  the  continuous-time 
autobicorrelation  or  autotricorrelation  that  is  sampled.  Generally,  it  is  the  original 
continuous-time  signal  that  is  sampled  and  used  to  calculate  the  discrete-time  autobicorrelation 
or  autotricorrelation.  In  this  case,  to  prevent  aliasing,  the  sampling  interval  for  the 
autobicorrelation  must  be  no  greater  than  two-thirds  the  interval  associated  with  critical 
sampling,  and  no  greater  than  one-half  for  the  autotricorrelation.  Numerical  calculations  of 
autocorrelation,  autobicorrelation,  and  autotricorrelation  zero-lag  values  corresponding  to 
spectral  area  and  volume  as  well  as  bispectral  contour  plots  for  two  model  bandlimited  energy 
signals  are  presented  as  demonstrations  of  these  conclusions.  Sampling  requirements  for 
higher-order  correlations  of  rectified  signals  are  also  discussed.  If  the  signal  is  to  be  rectified  for 
use  in  a  correlation  detector  or  time  delay  estimator,  then  the  sampling  rate  must  increase  to 
accommodate  the  higher  frequencies  that  usually  result  from  the  process  of  rectification. 
Spectral  masking  filters,  which  remove  aliasing  in  higher-order  correlations,  are  defined  and  a 
bispectral  filter  is  applied  to  some  representative  energy  signals.  Lag  domain  convolution  filters 
for  the  removal  of  aliasing  in  the  bicorrelation  and  tricorrelation  are  also  given.  These  filters 
assume  critical  sampling  or  finer  for  the  original  signal  to  remove  aliasing  effects  totally 
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*  _  Higher-order  correlations  and  spectra  may  be  used  to 

Wv  extract  information  about  signals,  which  the  autocorrela¬ 
tion  or  energy /power  spectrum  does  not  contain,  such  as 
phase,  deviations  from  Gaussianity,  and  the  presence  of 
nonlinearities  in  a  signal.  In  addition,  the  special  properties 
associated  with  higher-order  correlations  and  spectra  in- 
troduce  a  myriad  of  uses  for  them  in  the  areas  of  signal 
detection  (Nagata,  1978;  Hinich,  1982;  Dwyer,  1984,  1985; 
Brockett  etal.,  1987;  Hinich  et  al,,  1989;  G.  loup  etal., 
1989a,b;  Sullivan  and  Hinich,  1990;  Hinich  and  Wilson, 
1990;  Pike  etal.,  1991;  Pflug  etal.,  1992b),  time  delay 
estimation  (Nikias  and  Pan,  1988;  Pflug  etal.,  1990, 
1993),  and  phase  retrieval  (Bartelt  etal.,  1984;  Matsuoka 

*  —  and  Ulrych,  1984;  Nikias  and  Raghuveer,  1987),  which  in 
•  ^^^^tum  may  be  used  for  target  discrimination,  localization, 

; '  ^«^lassification,  etc.  For  most  practical  applications,  the  sig- 

"nal  to  be  studied  is  first  properly  digitized,  i.e.,  the  signal  is 
|«^enerally  bandlimited  according  to  the  sampling  rate  and 
then  sampled.  For  energy  signals,  direct  calculation  of 
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higher-order  correlations  using  a  digitized  signal  requires 
finer  sampling  relative  to  the  highest  frequency  present 
than  the  sampling  requirements  of  the  autocorrelation  to 
avoid  aliasing,  in  the  absence  of  a  masking  filter  (Pflug 
etal.,  1992a). 

In  Secs.  I-III,  definitions  and  Oescriptions  of  aliasing 
in  the  autocorrelation,  autobicorrelation,  and  autotricorre¬ 
lation  are  presented.  A  discussion  of  the  differences  in  sam¬ 
pling  requirements  for  (a)  discrete-time  correlations  from 
sampled  data  versus  (b)  those  obtained  by  sampling 
continuous-time  correlations  is  given  in  Sec.  IV.  In  Sec.  V, 
examples  of  critically  sampled  energy  signals  with  aliasing 
present  in  the  autobicorrelation  are  shown  and  tables  of 
spectral  area  and  volume,  as  defined  by  correlation  central 
ordinate  values,  are  given  to  demonstrate  the  sampling  re¬ 
quirements  for  second-,  third-,  and  fourth-order  correla¬ 
tions.  Higher-order  masking  filters  are  defined  in  Sec.  VI 
and  a  bispectral  masking  filter  is  applied  to  remove  aliasing 
in  the  autobicorrelation  of  the  signals  introduced  in  the 
previous  section,  which  for  this  application  are  critically 
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sampled.  Finally,  in  Sec.  VII,  the  effects  of  rectification  on 
signal  sampling  requirements  are  discussed. 


I.  ALIASING  IN  THE  AUTOCORRELATION 

Let  x(.t)  represent  a  real  continuous-time  function 
with  continuous-frequency  Fourier  transform  (Bracewell, 
1986) 


x(/)exp(  —i2irft)dt. 


The  inverse  transform  is  then 


x(t)=  X(f)e\p(i2Trft)df. 

«/  —  <x 

Existence  of  the  Fourier  transform  is  discussed  by  Papoulis 
(1962)  and  Bracewell  (1986).  If  x(f)  is  sampled  in  time 
with  sampling  interval  At  [resulting  in  the  discrete-time 
Fourier  transform  (DTFT)],  and  then  a  time  window  is 
applied  and  the  transform  is  sampled  in  frequency  with  N 
points,  the  result  is  the  discrete-time  Fourier  series 
(DTFS)  pair  given  by  Marple  (1987) 

.v-i 

X(k  Af)=At  2  x(nAt)cxp[—i2ir(n  At)(k  A/)/N] 

n  =  0 

(k=-N/2,...,-N/2-\), 
with  inverse 

/V/2-1 

x{nAt)=Af  X  X(k  Af)e\p[i2iT(nAt){kA/)/N] 

k=-N/2 


(n=0,...,n-l). 


and  A/=l/(,NAt).  Although  Marple  (1987)  distin¬ 
guishes  the  DTFS  from  the  DFT,  for  the  remainder  of  this 
paper  these  forms  will  be  referred  to  as  the  DFT  pair,  and 
we  use  t  in  place  of  n  At  and /in  place  of  A:  A/  in  discrete 
equations.  In  practice,  a  frequency  window  should  be  ap¬ 
plied  before  the  original  time  domain  signal  is  sampled.  If 
the  signal  is  bandlimited  in  frequency  and  the  continuous¬ 
time  function  is  oversampled,  i.e.,  /n=1/(2A/)  is  chosen 
such  that  /N>/r,  where  /n  represents  the  Nyquist  or 
foldover  frequency,  and  /,  is  the  top  or  highest  frequency 
present  in  the  signal,  then  the  centers  of  the  replications  of 
the  continuous  transform  are  spaced  greater  than  a  dis¬ 
tance  of  2/,  apart,  and  do  not  overlap.  In  this  case,  the 
principal  replication  may  be  isolated  using  a  rectangular 
frequency  window  defined  by  Bracewell  (1986), 


mf)= 


(1.  I/I<A 

lo,  |/|>/„ 


tion  is  not  an  accurate  representation  of  the  true 
continuous-time  transform,  that  is,  aliasing  has  occurred. 
In  other  words,  if  the  function  is  sampled  with  sampling 
interval  At>  At^,  where  At^  is  the  sampling  interval  nec¬ 
essary  for  critically  sampling  a  signal,  then  the  transform 
replications  overlap  in  the  principal  domain  and  the  trans¬ 
form  at  each  point  of  overlap  will  be  the  sum  of  the  una¬ 
liased  and  overlapping  transforms. 

To  avoid  confusion,  we  must  make  a  clear  distinction 
between  the  Nyquist  sampling  frequency,  /n,  and  the 
highest  or  top  frequency  present  in  a  bandlimited  signal, 
denoted  /,.  Whereas  /n  is  purely  a  characteristic  of  the 
sampling  interval,  /,  is  purely  a  function  of  the  individual 
signal.  For  an  essentially  naturally  bandlimited  signal, 
these  values  are  completely  independent  of  each  other. 
However,  for  a  noise  signal  with  theoretically  infinite  band¬ 
width,  /,  is  dependent  on  the  initial  sampling  rate,  al¬ 
though  the  signal  may  be  interpolated  (or  modified  with  a 
masking  filter  as  we  discuss  later)  to  avoid  aliasing  when 
using  the  digitized  signal  in  higher-order  correlation  cal¬ 
culations.  The  calculation  of  the  autocorrelation  of  a  sam¬ 
pled  signal  requires  only  the  bandwidth  of  the  original  sig¬ 
nal  to  avoid  aliasing.  Therefore,  if  a  discrete  energy  signal 
xU)  has  autocorrelation, 

1 

at{T)=  £  x{t)x(t-^T)At, 

1=0 

and  energy  spectrum  defined  by 

,v-i 

ESA/)=X{f)X*(f)=  I  a,(?)exp(-t77r//)Ar, 

f=0 

then  x(t)  must  be  sampled  with  sampling  interval 
A/<A/c=  1/(2/,)  to  avoid  aliasing  in  the  autocorrelation 
(Bracewell,  1986). 


II.  THE  BICORRELATION 
A.  Aliasing  in  the  bicorreiation 

The  autobicorrelation  and  corresponding  bispectrum 
of  random  signals  are  defined  in  Brillinger  and  Rosenblatt 
(1967a)  and  Hasselmann  et  al.  (1963).  The  autobicorre¬ 
lation  of  an  energy  signal  is  defined  as  a  natural  extension 
of  the  autocorrelation  of  an  energy  signal  (Bracewell, 
1986;  G.  loup  et  al.,  1989a;  Pflug  et  al,  1992a,b),  and  for 
continuous-time  x{t)  is  given  by 

a2(Ti,T2)=  j  X{t)x{t  +  Ti)x[t  +  T2)dt, 

and  for  discrete- time  x(r)  by 


for  which  ±  /,  define  the  minimum  passband  filter  limits. 
However,  the  limits  of  the  filter  can  be  any  frequency  up  to 
±(2/n— ft).  If  the  continuous-time  function  is  bandlim¬ 
ited  at  /,  with  the  frequency  content  at  ±  /,  finite,  and  the 
function  is  sampled  at  /n=//,  then  the  centers  of  the 
replications  in  the  frequency  domain  are  spaced  exactly  2 /, 
apart  (critical  sampling)  and  only  overlap  at  the  endpoints 
of  each  replication.  No  aliasing  occurs.  However,  if  the 
function  is  not  sufficiently  sampled,  the  principal  replica¬ 


/v-i 

02(.Ti,T2)=  X  X{t)x(t  +  T^)x{t  +  T2)At. 

t=0 

Its  two-dimensional  continuous-time  Fourier  transform, 
the  bispectrum,  is  defined  by  (Pflug  et  al,  1992a) 

B{f\,f2)=  \  r  a2{Tx,T2) 

J  —  00  •/  —  oc 

Xexp[-/2ir(Ti  f  \-^t2  f 2)\dt\  di-j. 
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and  in  discrete  form, 

JV-I  V-l 

B(f\,h)=  X  X  02(^1  .r2)exp[-/2ir{T,/, 

ri=0  r2=0 

+  T2/2)]At,  AT2. 

The  bispectrum  of  the  discrete-time  autobicorrelation  is 
replicated  with  period  2/^  in  the  /1-/2  plane;  however,  the 
principal  domain  or  principal  replication  may  be  defined 
by  restricting  /i  and  sucn  that— /n</i</n  and 
— /n</2</n- 

It  has  been  shown  that  the  bispectrum  of  a  real  energy 
transient  may  be  expressed  as  a  product  of  one-dimensional 
Fourier  transforms  (Pflug  et  ai,  1992a),  that  is, 

fx,fi)=Xi  f^)X{  /,  +  /2) 


or 


5(  /,,/2)=WiW/2W3). 

where  f i=  —  f \—  f i-  The  three-index  notation  was  origi¬ 
nally  introduced  by  Brillinger  and  Rosenblatt  ( 1967b)  and 
discussed  in  Rosenblatt  (1985)  for  stationary  random  pro¬ 
cesses.  However,  it  applies  to  finite  length  energy  transients 
as  well  (G.  loup  et  ai,  1989a,b,  1990;  Pflug  etal, 
1992a,b).  A  derivation  of  this  expression  for  energy  signals 
is  given  by  Pflug  etal  (1992a).  The  total  unaliased  do¬ 
main  (TUD)  of  the  autobispectrum  of  a  sampled  signal  is 
the  result  of  the  overlap  of  the  domains  of  the  two  factors 
X{f^)X{f2)  and  (G.  loup  etal,  1989a, 

1990;  Pflug  et  ai,  1989,  1992a).  When  the  autobicorrela¬ 
tion  is  calculated  from  discrete-time  data,  the  factors  in  the 
above  equation  for  Bif^,/^)  must  be  replicated  infinitely 
with  period  2/^  in  the  bifrequency  plane,  i.e., 

5(/„/2)  =  f  i  Jr(/,+2A:,/N)) 

V*i  =  -«  / 

xf  X  X(  fi+lki  f^) 

\  *2=-« 

xf  X  X*(  fi+f2  +  2kj  ffi) 

\  *3= -00 

This  result  is  derived  in  Appendix  A  for  the  bispectrum  of 
the  bicorrelation  calculated  from  discrete-time  energy  sig¬ 
nals.  The  corresponding  result  for  discrete-time  stationary 
stochastic  signals  is  derived  in  Appendix  B. 

The  total  domain  of  the  bispectrum  results  from  the 
overlap  of  each  individual  term  in  each  summation  of  the 
above  equation,  as  shown  for  a  small  part  of  the  total 
bifrequency  domain  by  the  shaded  regions  in  Fig.  1  for  a 
signal  with  bottom  frequency  /*=0  and  /n= /»•  The  anal¬ 
ysis  of  the  overlap  is  facilitated  by  considering  the  domains 
of  the  replications  of  X(f^)X(f2)  separately  from  the  do¬ 
mains  of  the  replications  of  X{—f\—  f2)  (G.  loup  etai, 
1989a;  Pflug  et  ai,  1992a).  In  Fig.  1  the  factors  that  con¬ 
tribute  to  the  base  replication  are  shaded,  with  the  princi¬ 
pal  replication  of  the  TUD  outlined.  The  area  in  the  prin¬ 
cipal  replication  outside  the  TUD  of  the  bispectrum  in  Fig. 


FIG.  1.  Domain  of  the  bispectral  factors  which  overlap  to  form  the  TUD 
and  TAD  of  the  principal  replication  of  the  bispectrum  when  a  signal  is 
sampled  such  that  /n=/,  • 


1  is  the  total  aliased  domain  (TAD).  This  aliasing  in  the 
bispectrum  is  a  result  of  neighboring  diagonal  factors 
.T(—  /|— /2±2/n)  overlapping  the  principal  replication, 
similar  to  aliasing  in  the  ordinary  spectrum. 

Figure  1  shows  aliasing  in  the  bispectrum  as  a  result  of 
critically  sampling  the  original  x(/).  While  the  type  of 
aliasing  caused  by  overlap  of  the  factor  X(fx)X{f2)  and 
the  replications  of  the  factor  A’(  —  f\—f 2)  to  the  immedi¬ 
ate  left  and  right  of  the  principal  diagonal  replication  is 
removable  via  a  masking  filter  (which  is  discussed  in  more 
detail  later),  sampling  below  the  critical  rate  results  in 
additional  overlap  of  the  factor  X(fx)X{f2)  with 
-T(/i±2/n)A'(/2±2/n)  as  well  as  the  diagonal  factors, 
causing  aliasing  that  is  not  removable.  To  avoid  aliasing  of 
either  kind  in  the  bispectrum  of  a  sampled  signal,  the  cen¬ 
ters  of  the  replications  must  be  a  distance  at  least  3  f,  apart. 
For  the  energy  spectrum,  if  /,  is  the  top  frequency  present 
in  a  signal,  then  A/c<l/(2/,);  but  for  the  bispectrum  the 
sampling  interval  Atg  must  be  (Pflug  et  ai,  1992a) 

Af2,<l/(3/,)  =  (2/3)Atc. 

If  the  signal  is  sampled  such  that  /,</n<3//2,  then  the 
aliasing  can  be  removed  by  using  a  masking  filter,  as  dis¬ 
cussed  in  Sec.  VI. 

B.  Symmetries  in  the  bispectrum 

Permuting  arguments  in  the  expression 
B{f\,f2)=X{f\)X{f2)X{—f^-f2)  leaves  the  value  of 
the  bispectrum  unchanged.  This  results  in  six  lints  of  sym¬ 
metry  which  intersect  the  principal  replication  of  the 
bispectrum  (Brillinger  and  Rosenblatt,  1967b): 

/i=/2>  2/i-f/2=0,  /i-f- 2/2=0, 

/,=0,  /2  =  0,  /,  =  -/2. 

The  last  three  symmetry  lines  involve  a  complex  conjugate 
operation.  The  six  symmetry  lines  divide  the  TUD  into  12 
triangular  domains,  each  containing  the  same  bispectral 
values.  From  the  more  general  expression  given  above  for 
the  bispectrum  of  a  bicorrelation  calculated  from  sampled 
data,  we  find  four  more  symmetry  lines  which  intersect  the 
principal  replication  of  the  bispectnim,  that  is.  for 
^]=^2=0  and  ^3=  ^  1,  tiiC  aliased  parts  of  the  bispectrum 
are  divided  by  the  lines 
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FIG.  2.  Symmetry  lines  in  the  TUD  and  TAD  of  the  bispectrum  calcu¬ 
lated  from  a  discrete  signal  sampled  such  that  /n  =  /,  •  Symmetry  lines  in 
the  TUD  apply  to  the  bispectrum  of  a  continuous-time  autobicorrelation 
as  well. 

2/i+/2=  ±2/n 
and 

f\  +  2/2  =  ±  2/n  • 

Like  the  TUD,  the  total  aliased  domain  is  divided  into  12 
regions.  The  axes  and  the  dashed  lines  shown  in  Fig.  2  are 
the  lines  of  symmetry  present  in  the  bispectrum  when  a 
signal  is  sampled  such  that  /n= /,  (/a=0  in  this  figure). 
The  entire  bispectrum  may  be  mapped  from  values  in  the 
principal  unaliased  polygon  (PUP)  and  the  principal 
aliased  polygon  (PAP)  using  these  symmetries  (Nikias 
and  Raghuveer,  1987;  Hinich  and  Wolinsky,  1988;  Pflug 
etaL,  1992a). 

The  symmetries  in  the  bispectrum  are  a  function  of  the 
sampling,  i.e.,  a  function  of  /n  ,  regardless  of  the  f,  for  the 
individual  signal.  The  bispectral  symmetries  are  shown  in 
Fig.  3  for  a  signal  with  0</*<  f,/l  and  sampled  with 
chosen  to  be  greater  than  /, .  If  the  signal  is  sampled  such 
that  /n=/;.  then  the  bispectral  domain  is  slightly  differ¬ 
ent,  as  shown  in  Fig.  4.  In  this  figure,  PAP  is  a  four-sided 
polygon,  not  a  triangle  as  shown  in  Fig.  3  with  /n>/»- 


FIG.  3.  Domain  of  the  bispectrum  of  a  discrete  signal  with  0<  /j<;//2, 
sampled  such  that  /n  >  /, . 


FIG.  4.  Domain  of  the  bispectrum  of  a  discrete  signal  with  0<  //,</  '2,  • 

sampled  such  that  fs  =  f, 

The  TUD  disappears  for  /*>  /,/2,  in  general.  If  is 
chosen  such  that  if i/2<ff^<.3f,/2  and  fi,>  f,/2,  then  the 
TUD  of  the  bispectrum  disappears  while  the  TAD  is  still 
present  (Pflug  e/ a/.,  1992a). 

III.  THE  TRICORRELATION 

A.  Aliasing  in  the  tricorretation 

To  discuss  aliasing  in  the  tricorrelation,  we  follow  rea¬ 
soning  similar  to  that  given  for  the  bicorrelation.  The 
continuous-time  autotricorrelation  03  (tj  ,72,73)  and  its 
three-dimensional  Fourier  transform  are  given  by  (Brill- 
inger  and  Rosenblatt,  1967a;  G.  loup  et  ai,  1989a) 

“3(71,72,73)=  I  x(f)x(t-(-7|)x(f-t-T2)x(r-t-73)df 

«/  —  cc 

and 

n/i,/2./3)=  r  r  f  “3(71,72,73) 

J  —  cc  J  —  oc  J  —  cc 

Xexp[  -|■277(7,/,-^-T2/2 
+  rifi)]dTx  dTidTi, 

and  the  discrete-time  autotricorrelation  and  trispectrum 

are  given  by  ^ 

A- I 

“3(71,72,73)=  X  x(r)x(f-f T,)x(r-l-r2)x(f-|-T3)Af 

/=0  V 

and 

A-_i  N-\  yv-i 
X  2  X  “3(71,72,73) 

r|=0  r2  — 0  T^  =  0 

X  exp  [  -  i2ir{  7,  /,  -I-  73  fi 

+  73/3)^71  A72Ar3.  ^ 

The  trispectrum  of  a  real  finite  energy  transient  may  equiv¬ 
alently  be  expressed  as  (Brillinger  and  Rosenblatt,  1967b; 

Dalle  Molle  and  Hinich,  1989;  Pflug,  1990;  Pflug  et  al., 

1992a) 
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T{  /,.  /2.  h)=X{  fOX{  h)X(  f,)X*{  /,+/2 


T{  fuf2,f,)=X(  fOX{  f^)X{  h)X{  /4). 

where  f —f  \— f 2— f  i  -  For  an  autotricorrelation  calcu¬ 
lated  from  a  sampled  signal,  each  factor  is  replicated  infi¬ 
nitely  in  the  trifrequency  plane,  i.e., 

F(/„/2./3) 

-I  i  Jr(/,+2A:,/N)W  i  X(f2 

=  —  X  /  \  ^2~  ~  ^ 

+  2^2/N)j^^  S  ■^(  /3T2/C3 /n)  j 


x*{  /l -H /2  4-/3  + 2^4 /n) 


The  TUD  of  the  trispectrum  calculated  from  a  sampled 
signal  is  defined  by  the  overlap  of  the  nonzero  domain  of 
the  two  factors  X{fi)X{f2)X(fi)  and  -1-/2 -1-/3) 

in  the  principal  replication  as  discussed  by  Pflug  et  al. 
( 1992a).  The  remaining  area  in  the  principal  replication  is 
a  result  of  the  overlap  of  the  factor 

x{  A)Xi/2)Xif2) 

from  the  principal  replication  and 

X(-A-/2-A-2ka^) 

with  k4=  ±  1,  and  comprises  the  TAD  of  the  trispectrum. 

To  avoid  the  kind  of  aliasing  caused  by  sampling  a 
signal  with  /n=//.  the  centers  of  the  trispectral  replica¬ 
tions  must  be  a  minimum  distance  of  4 /,  apart.  Hence,  the 
sampling  interval  Atj-  must  be  (Pflug  etal.,  1992a) 

Atr<  1/(4/,)  =  (1/2)  Ate- 

If  the  signal  xit)  is  sampled  such  that  /,</n<2  /,,  then 
the  aliasing  in  the  trispectrum  can  be  removed  with  a 
masking  filter  as  discussed  below. 

B.  Symmetries  in  the  trispectrum 

Symmetry  planes  in  the  trispectrum  are  derived  by 
permuting  the  arguments  in  the  above  expression  for  the 
trispectrum  of  a  tricorrelation  calculated  from  sampled 
data.  When  k^=k2  =  k2  =  k^=0,  the  symmetry  planes  are 
(Brillinger  and  Rosenblatt,  1967b;  Dalle  Molle  and  Hin- 
ich,  1989;  Pflug,  1990;  Pflug  et  al.,  1992a) 

f\=A>  A  =  fi’  /l=/3.  2/i4-/2  4-/3  =  0, 

/|+2/2  +  /3  =  0>  /i+/2  +  2/3  =  0, 

A=0,  /2  =  0,  /3  =  0. 

These  planes  divide  the  TUD  of  the  trispectrum  into  two 
types  of  unaliased  principal  polyhedral  domains,  PUPl 
and  PUP2,  each  replicated  48  times.  PUPl,  as  defined  by 
Pflug  etal  (1992a),  disappears  for  /*> /,/3.  The  TAD 


also  contains  two  types  of  principal  aliased  polyhedra, 
PAPl  and  PAP2,  a  result  of  the  intersection  of  the  sym¬ 
metry  planes 

2/i+/2  +  /3=  ±2/n, 

/i  +  2/2+ f  I  =  ±  2/n  , 

2/3  =  ±  2/n  , 

with  the  TAD.  Both  PAPl  and  PAP2  are  each  replicated 
48  times  in  the  TAD.  Detailed  discussion  and  figures  de¬ 
picting  the  TUD,  TAD,  PUP,  and  PAP  of  the  trispectrum 
can  be  found  in  Pflug  et  al  ( 1992a). 

IV.  SAMPLING  THE  DATA  VERSUS  SAMPLING  THE 
CORRELATION 

The  characterization  of  aliasing  for  functionals  such  as 
the  correlation  contains  a  subtlety  that  is  often  unrecog¬ 
nized.  That  is,  the  nature  of  the  aliasing  and  even  the  re¬ 
quirements  to  prevent  it  can  depend  on  whether  the  cor¬ 
relation  is  calculated  from  data  that  are  already  sampled  or 
the  continuous-time  correlation  itself  is  sampled.  Follow¬ 
ing  the  reasoning  of  Kanasewich  ( 1973)  the  expression  for 
the  transform  of  the  autocorrelation,  which  describes  the 
effects  of  aliasing  for  a  sampled  continuous-time  autocor¬ 
relation,  is 

00 

ES,{f)=  2  \X(f-lkf^)\\ 

k~  —  ao 

where 

oc 

X(f)=  2  X{f-2kf^) 

k=  —  cc 

is  the  aliased  Fourier  transform.  The  tilde  is  used  to  indi¬ 
cate  aliased  functions  and  transforms.  This  result  does  in¬ 
deed  describe  the  aliasing  if  the  continuous-time  autocor¬ 
relation  itself  is  sampled.  If,  however,  the  discrete-time 
autocorrelation  is  calculated  from  data  x(t),  which  have 
already  been  sampled  such  that  x(t)  contains  aliasing,  then 
the  appropriate  transform  domain  expression  correspond¬ 
ing  to  the  aliasing  of  the  autocorrelation  is 

_ ^00  2 

£SA/)=  2  X(f-lkf2i)  . 

k=  —  OD 

The  difference  is  that  the  unaliased  transform  in  the 
first  equation  is  replicated  after  squaring  whereas  the  una¬ 
liased  X{f)  in  the  second  equation  is  replicated  before 
squaring,  giving  rise  to  interference  due  to  phase.  The  re¬ 
quirement  to  prevent  aliasing,  namely  that  /N>/n  with  /, 
the  highest  frequency  present  in  the  unsampled  data  x(t), 
is  the  same  for  both  cases,  however. 

The  characterization  of  aliasing  for  higher-order  cor¬ 
relations  contains  an  additional  striking  feature.  The  un¬ 
sampled  continuous-time  autobicorrelation  transform, 
X(/i)X(,A)X*{f \+A)y  has  a  nonzero  domain  which 
can  be  determined  from  the  overlap  of  the  nonzero  do¬ 
mains  of  the  individual  factors.  This  method  has  been  used 
previously  by  G.  loup  et  al  ( 1989a,b  and  1990),  and  Pflug 
etal  (1992a).  The  nonzero  domains  of  continuous-time 
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FIG.  5.  TUD  of  the  bispectrum  of  a  continuous-time  autobicorrelation. 

transforms  1 )  and  A'(/2)  are  the  infinite  strips  defined 
by  and  — /,</2</,.  respectively,  and  the 

nonzero  domain  of  A'*(  /i  -I-/2)  is  the  infinite  strip  defined 
by  —  /,</i  +  /2<//-  These  are  shown  in  Fig.  5.  The  over¬ 
lap  of  these  three  nonzero  domains  is  the  heavily  shaded 
region  in  Fig.  5,  which  corresponds  to  the  TUD  pictured  in 
Figs.  1  and  2.  All  the  nonzero  autobispectrum  (the  trans¬ 
form  of  the  continuous-time  autobicorrelation)  is  con¬ 
tained  in  a  square  that  is  2/,  on  a  side,  so  the  sampling 
required  to  prevent  aliasing  if  the  continuous-time  autobi¬ 
correlation  itself  is  sampled  is  /n>/i.  the  same  as  that 
required  for  the  function  and  the  autocorrelation.  This  re¬ 
sult  has  been  derived  independently  by  Nielsen  (1992). 
Pflug  et  al.  (1992a)  have  shown  that  the  sampling  of  the 
data  needed  to  prevent  aliasing  for  the  discrete-time  auto¬ 
bicorrelation  calculated  from  sampled  data  is  /n>3//2. 
An  extension  of  this  argument  can  be  made  for  sampling 
the  continuous-time  autotricorrelation  directly,  for  which 
the  requirement  to  prevent  aliasing  is  also  /N>/f  If  sam¬ 
pled  data  are  to  be  used,  the  requirement  is  Note 

that  there  will  be  no  TAD  region  if  a  continuous-time 
higher-order  correlation  is  sampled  directly. 

It  is  important  to  consider  the  applicability  of  these 
two  types  of  models.  Most  calculations  of  the  correlation 
involve  the  use  of  sampled  data,  so  the  models  for  aliasing 
in  the  correlations  due  to  calculations  from  sampled  data 
are  most  often  applicable.  The  sampling  of  the  correlation 
directly  occurs  much  less  frequently  and  principally  when 
the  correlation  itself  is  the  result  of  a  physical  measure¬ 
ment.  For  these  cases  the  model  for  directly  sampling  the 
continuous-time  correlation  applies. 

The  obvious  difference  between  sampling  the 
continuous-time  correlation  and  calculating  the  bicorrela¬ 
tion  from  sampled  data  is  the  presence  in  the  latter  situa¬ 
tion  of  what  is  referred  to  in  this  paper  as  the  TAD.  Ali¬ 
asing  in  the  bicorrelation  and  tricorrelation  due  to 
sampling  such  that  /,</n<3//2  and  /,</n<2/,,  respec¬ 
tively,  is  removable  via  a  masking  filter,  while  sampling 
such  that  f^<f,  results  in  aliasing  that  is  not  generally 


removable.  One  might  choose  not  to  refer  to  the  TAD  of 
the  bispectrum  and  trispectrum  as  the  aliased  regions  since 
they  are  completely  removable  and  only  exist  when  the 
higher-order  correlation  is  calculated  from  data  that  are 
undersampled  for  this  purpose.  Hinich  et  al.  ( Hinich  and 
Wolinsky,  1988;  Hinich  et  al,  1989;  Hinich  and  Wilson, 

1990)  do  not  refer  to  the  contents  of  this  region  as  aliasing, 
but  instead  call  the  principal  domain  of  this  region  the  OT. 
Hasselmann  et  al.  (1963)  use  the  term  aliasing  to  refer  to 
this  region  of  the  bispectrum,  and  Brillinger  and  Rosen¬ 
blatt  (1967b)  discuss  these  regions  of  the  bispectrum  and 
trispectrum  in  the  section  of  their  paper  on  aliasing.  Our 
reasons  for  using  the  name  aliasing  for  the  TAD  are  given  ^ 
in  the  following  discussion.  •  ^ 

Aliasing  is  the  lack  of  sufficient  sampling  to  represent 
correctly  the  underlying  continuous-time  function  caused  ^ 
by  higher  frequency  content  folded  over  onto,  and  mas¬ 
querading  as,  lower  frequency  content.  Sampling  data  such 
that  /,</n  <  3 /,/2  for  a  bicorrelation  calculation  or  such 
that  /,</n  <2/,  for  a  tricorrelation  calculation  causes  the 
resulting  discrete-time  correlations  to  be  inaccurate  repre¬ 
sentations  of  the  continuous-time  correlations,  due  to  dis¬ 
tortion  of  the  Fourier  transforms  in  their  two-  or  three- 
dimensional  principal  domains,  which  result  from  the 
replication  of  the  transforms  produced  by  sampling  the 
continuous-time  data,  as  is  the  case  for  ordinary  aliasing. 
Moreover,  the  filtering,  which  is  normally  done  to  recon¬ 
struct  the  continuous-time  two-  or  three-dimensional  func¬ 
tion,  i.e.,  excluding  all  frequencies  outside  the  principal 
domain  or  replication,  does  not  remove  the  TAD.  It  is 
possible  to  argue  that  a  special  filter  should  be  used  to 
reconstruct  the  higher-order  correlation,  indeed,  the  filter 
that  we  discuss,  but  this  is  not  in  keeping  with  general 
reconstruction  for  arbitrary  two-  or  three-variable  func¬ 
tions.  An  additional  important  feature  is  that  while  the 
autocorrelation  calculated  from  data  sampled  such  that 
/n=//  is  correct  at  the  sample  points,  the  higher-order 
autocorrelations  calculated  from  the  same  data  are  incor¬ 
rect  at  the  sample  points.  For  these  important  reasons,  the 
regions  of  the  bispectrum  and  trispectrum  defined  by  sam¬ 
pling  such  that  /,</n<3//2  and  /,</n<2/,,  respec¬ 
tively,  are  referred  to  as  aliased  domains. 

V.  EXAMPLES  OF  ALIASING  IN  THE  CORRELATIONS 
OF  SAMPLED  SIGNALS 

Two  model  energy  signals  are  used  to  show  the  occur-  4 
rence  of  aliasing  in  higher-order  correlations  and  spectra 
due  to  sampling  a  signal  at  /n=/c  The  first  is  a  model  of 
an  amplitude  and  frequency  modulated  signal  for  a  20-Hz 
finback  whale  signal  fitted  by  J.  loup  et  al.  ( 1988;  G.  loup 
et  aL,  1989a)  and  modeled  after  a  real  signal  shown  in 
Watkins  et  al.  ( 1987),  with  frequencies  equal  to  and  above 
32-Hz  set  to  zero  for  simplicity.  The  signal  and  its  Fourier 
magnitude  spectrum  are  shown  in  Fig.  6.  Notice  that  am¬ 
plitudes  for  frequencies  below  about  10  Hz  are  quite  small.  * 
Hence  the  signal  has  a  narrow-band  character  with  /,  not 
much  larger  than  2/*,  if  /*  is  considered  to  be  near  10  Hz. 

Thus  we  expect  the  volume  under  its  bispectrum  to  be 
small  (Pflug  et  al.,  1992a).  Applying  the  conditions  de- 
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FIG.  6.  Simulated  whale  transient  and  its  Fourier  magnitude  spectrum. 

scribed  in  earlier  sections,  the  signal  should  be  sampled 
with  /n=32  Hz  to  avoid  aliasing  in  the  autocorrelation, 
with  /n  =  48  Hz  to  avoid  aliasing  in  the  autobicorrelation, 
and  with  /n  =  64  Hz  to  avoid  aliasing  in  the  autotricorre- 
lation.  For  the  whale  signal  of  1-s  duration,  this  corre¬ 
sponds  to  64  sample  points  to  avoid  aliasing  in  the  auto¬ 
correlation,  96  sample  points  for  the  autobicorrelation,  and 
128  points  for  the  autotricorrelation. 

Ideally,  if  there  is  no  aliasing  in  a  correlation  calcu¬ 
lated  from  sampled  data,  finer  sampling  in  a  fixed  time 
period  will  not  change  the  zero-lag  value,  i.e..  the  area 
beneath  the  energy /power  spectrum  (Kanasewich,  1973). 
Other  lag  values  could  be  used  to  demonstrate  these  prop¬ 
erties.  The  whale  signal  is  sampled  with  various  numbers  of 
points  and  the  area  beneath  the  energy  spectrum  and  vol¬ 
umes  beneath  the  bispectrum  and  trispectrum  are  calcu¬ 
lated  by  zero-lag  correlation  value.  Results  to  five  signifi¬ 
cant  digits  are  shown  in  Table  I.  The  zero  lag  value  of  the 
rectified  autobicorrelation  is  also  given,  and  will  be  dis¬ 


cussed  later.  While  the  autocorrelation  requires  only  ap¬ 
proximately  64  sample  points,  as  predicted,  for  the  effects 
of  aliasing  to  become  negligible  to  five  significant  digits,  the 
autobicorrelation  and  autotricorrelation  each  require  more 
sample  points.  The  autobicorrelation  requires  the  predicted 
96  sample  points  before  the  aliasing  disappears,  and  the 
hypervolume  beneath  the  autotricorrelation  is  constant  to 
five  significant  digits  if  a  minimum  of  104  points  is  used, 
somewhat  fewer  than  predicted. 

A  color  contour  plot  of  the  magnitude  bispectrum  of 
the  whale  signal  sampled  with  64  points  is  shown  in  Fig.  7 
with  the  bispectral  symmetry  lines  distinguishing  the  una¬ 
liased  and  aliased  regions  (compare  to  Fig.  2).  As  we  pre¬ 
dict,  the  bispectrum  contains  aliasing.  However,  the  energy 
spectrum  with  a  nonzero  amplitude  at  a  maximum  of  32 
Hz  (see  Fig.  6  for  the  Fourier  magnitude  spectrum)  con¬ 
tains  no  aliasing  at  this  rate  of  sampling.  If  128  sample 
points  are  used,  the  aliasing  in  the  bispectrum  disappears, 
as  shown  in  Fig.  8,  which  has  a  different  normalization 
from  that  used  for  Fig.  7. 

In  contrast  to  the  narrow-band  whale  signal,  the  sec¬ 
ond  example  is  a  narrow  pulse  (narrow  in  time,  broadband 
in  frequency)  shown  with  its  Fourier  magnitude  in  Fig.  9. 
For  this  signal  of  1  -s  duration,  we  expect  that  aliasing  will 
not  exist  in  the  autocorrelation  if  the  signal  is  sampled  such 
that  /n= /j=256  Hz,  or  with  512  sample  points.  To  avoid 
aliasing  in  the  autobicorrelation,  the  signal  must  be  sam¬ 
pled  with  at  least  768  points,  and  to  avoid  aliasing  in  the 
autotricorrehtion,  the  signal  must  be  sampled  with  at  least 
1024  points.  We  see  in  Table  II  that  the  zero  lag  value  of 
the  autocorrelation  is  constant  to  five  significant  digits  if  a 
minimum  of  512  sample  points  is  used.  The  autobicorrela¬ 
tion  central  ordinate  is  constant  to  five  significant  digits 
using  a  minimum  of  704  points,  and  the  autotricorrelation 
central  ordinate  is  constant  using  a  minimum  of  832 
points.  Using  the  predicted  768  points  or  more  for  the 


TABLE  1.  TimeKiomain  calculations  of  correlation  central  ordinate  values  of  the  whale  transient  with  various  sampling  rates. 


Number  of  points 

Autocorrelation 
(amp  units)^  s 

Whale  transient 
Autobicorrelation 
(amp  units)^  •  s 

Autotricorrelation 
(amp  units)*  ■  s 

Rectified  autobicorrelation 
(amp  units)'  ■  s 

32 

2.2296  xlO'^ 

3.4944x10-’ 

1.0268  XlO  ’ 

4.6164  xlO  ' 

48 

0.27676 

3.3260X10-^ 

0.17467 

0.21102 

64 

0.27708 

3.0329X10-’ 

0.16384 

0.20488 

72 

0.27708 

6.1817x10-“ 

0.16094 

0.20364 

80 

0.27708 

6.8907X10-“ 

0.16549 

0.20585 

88 

0.27708 

6.8659X  10-“ 

0.17463 

0.21190 

96 

0.27708 

6.8661X10  “ 

0.17413 

0.21093 

104 

0.27708 

6.6861X10-* 

0.17414 

0.21051 

112 

0.27708 

6.8661X10  “ 

0.17414 

0.21129 

120 

0.27708 

6.8661  X 10  “ 

0.17414 

0.21211 

128 

0.27708 

6.8661X10  “ 

0.17414 

0.21161 

136 

0.27708 

6.8661  X 10  “ 

0.17414 

0.21136 

200 

0.27708 

6.8661X10  “ 

0.17414 

0.21143 

300 

0.27708 

6.8661X10  “ 

0.17414 

0.21141 

400 

0.27708 

6.8661x10-* 

0.17414 

0.21142 

500 

0.27708 

6.8661x10  “ 

0.17414 

0.21142 

600 

0.27708 

6.8661x10  “ 

0.17414 

0.21142 
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FIG.  7.  Contour  plot  of  the  magnitude  bispectrum  of  the  whale  transient  FIG.  8.  Contour  plot  of  magnitude  bispectrum  of  the  whale  transient 

sampled  with  64  points  so  that  aliasing  is  present.  The  amplitudes  have  sampled  with  128  points  so  that  no  aliasing  is  present.  The  amplitudes 

been  normalized  to  unit  height.  have  been  normalized  to  unit  height. 
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FIG.  9.  Narrow  pulse  and  its  Fourier  magnitude  spectrum. 
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FIG.  10.  Contour  plot  of  the  first  quadrant  of  the  magnitude  bispectrum  FIG.  It.  Contour  plot  of  the  first  quadrant  of  the  magnitude  bispectrum 

of  the  narrow  pulse  sampled  with  S12  points  so  that  aliasing  is  present.  of  the  narrow  pulse  sampled  with  1024  points  so  that  no  aliasing  is 

The  amplitudes  have  been  normalized  to  unit  height.  present.  The  amplitudes  have  been  normalized  to  unit  height 


TABLE  II.  Time-domain  calculations  of  correlation  central  ordinate  values  of  the  narrow  pulse  transient  with  various  sampling  rates. 


Number  of  points 

Autocorrelation 

(XlO-’) 

(amp  units •  s 

Narrow  pulse 
Autobicorrelation 
(xio-’) 

(amp  units)^-s 

Autotricorrelation 

(xio-^) 

(amp  un-*s)*  •  s 

Rectified  autobicorrelation 
(XlO-^) 

(amp  units ■  s 

256 

2.3263 

1.7867 

1.3766 

1.7868 

320 

2.5966 

2.2879 

2.0588 

2.2892 

384 

2.6770 

2.4089 

2.3176 

2.4170 

448 

2.6875 

2.3150 

2.2257 

2.3308 

512 

2.6880 

2.2331 

2.0731 

2.2441 

576 

2.6880 

2.2055 

1.9787 

2.2100 

640 

2.6880 

2.2019 

1.9426 

2.2093 

704 

2.6880 

2.2017 

1.9324 

2.2188 

768 

2.6880 

2.2017 

1.9305 

2.2108 

832 

2.6880 

2.2017 

1.9303 

2.2094 

896 

2.6880 

2.2017 

1.9303 

2.2094 

960 

2.6880 

2.2017 

1.9303 

2.2103 

1024 

2.6880 

2.2017 

1.9303 

2.2104 

2000 

2.6880 

2.2017 

1.9303 

2.2101 

3000 

2.6880 

2.2017 

1.9303 

2.2100 

4000 

2.6880 

2.2017 

1.9303 

2.2100 

5000 

2.6880 

2.2017 

1.9303 

2.2100 

6000 

2.6880 

2.2017 

1.9303 

2.2100 

7000 

2.6880 

2.2017 

1.9303 

2.2100 

8000 

2.6880 

2.2017 

1.9303 

2.2100 
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autobicorrelation  and  a  double  precision  calculation,  the 
central  ordinate  value  is  constant  to  ten  significant  digits, 
and  the  autotricorrelation  central  ordinate  is  constant  to  14 
significant  digits  using  the  predicted  1024  points  or  more. 
In  each  case,  the  central  value  is  constant  to  ten  significant 
digits  using  no  more  than  the  predicted  number  of  sample 
points.  Comparing  the  first  quadrant  of  the  magnitude  bis¬ 
pectrum  of  the  narrow  pulse  sampled  with  512  points, 
shown  in  Fig.  10,  and  the  Fourier  magnitude  spectrum  of 
the  original  signal  shown  in  Fig.  9,  we  see  that  aliasing  is 
present  in  the  bispectrum  if  only  512  points  are  used  even 
though  the  energy  spectrum  contains  no  aliasing  if  sampled 
such  that  /n>256  Hz.  If  1024  points  are  used,  the  bispec¬ 
trum  no  longer  shows  any  aliasing,  as  shown  in  Fig.  11. 


VI.  HIGHER-ORDER  MASKING  FILTERS 


Aliasing  in  higher-order  correlations  due  to  sampling 
such  that  /n> /,  may  be  removed  using  one  of  two  tech¬ 
niques.  The  first  is  to  operate  on  the  original  discrete-time 
signal.  The  discrete-time  signal  may  be  interpolated  to  in¬ 
clude  a  sufficient  number  of  points  within  the  signal  dura¬ 
tion  to  avoid  aliasing  in  the  autobicorrelation  or  autotri¬ 
correlation.  This  is  most  often  done  by  zero-padding  its 
Fourier  transform. 

The  second  technique  involves  operating  on  the  fre¬ 
quency  domain  autobispectrum  or  autotrispectrum  by  ap¬ 
plying  a  higher-order  masking  filter  or  on  the  aliased  au¬ 
tobicorrelation  or  autotricorrelation  with  a  convolutional 
masking  filter.  A  bispectral  masking  filter  in  the  frequency 
domain  is  defined  by  (Pflug  et  ai,  1992a) 


1, 

0, 


1/i+/2|</m 


or  in  the  lag  domain, 


X 

h{T^,T2)=^f]  'L  sinc(2/^)5<nc[2/,(/-|-T,)] 

l~=  —  X 


Xsinc[2/,(r-f  rjllA/, 

where  the  time  domain  filter  h[T ^,72)  is  the  autobicorrela¬ 
tion  of  the  sine  function.  If  instead  of  the  continuous- 
frequency  filter,  the  time  domain  equivalent  of  a  discrete- 
frequency  filter  is  sought,  the  above  sum  becomes  a  finite 
sum  over  digital  sines  (Marple,  1987).  If  /n> /,,  then  any 
value  between  /,  and  may  be  used  as  the  bandlimit  in 
the  filter.  The  unaliased  bispectrum  is 

\,f 2)  =  Bif  \,f 2)H(f  2)  and  the  corresponding  un¬ 
aliased  autobicorrelation  02(^1  >’’2)  =  02(^1  >^2)*^(’'i»^2)- 
Frequency-domain  multiplications  are  more  commonly 
used  since  they  are  faster  than  time-domain  convolutions. 
Indeed,  convolving  the  two-dimensional  time  domain  filter 
/7(T|,T2)  with  the  aliased  autobicorrelation  could  be  pro¬ 
hibitive. 

A  bispectral  filter  //(/i,/2)  was  applied  to  the  whale 
transient  sampled  with  64  points,  a  number  just  sufficient 
to  prevent  aliasing  in  the  autocorrelation.  As  shown  in 
Table  I,  the  volume  beneath  the  bispectrum  is 
3.0329X  10“^  (amp  units)^  •  s  without  the  filter,  a  341.7% 


error  from  the  expected  unaliased  value  of  6.8661  X  10  * 
(amp  units)^  s.  However,  the  volume  is  reduced  to 
6.6208x10  (amp  units)^  s  when  the  filter  is  applied, 
resulting  in  a  3.7%  error.  For  the  narrow  pulse  sampled 
with  512  points,  a  masking  filter  reduced  the  aliased  vol¬ 
ume  from  2.2331  X  10~^  (amp  units)'  s,  a  1.4%  error,  to 
2.2017 X  10“'  (amp  units)'  s.  which  has  zero  error  up  to 
five  significant  digits  (see  Table  II). 

A  trispectral  masking  filter  can  be  written  as 


J  >  I  /i  +  /2  +  /a  I  . 

0.  I  /l  +  /2  +  /.^  I  >  /n 


or 


A(T,,T2,rj)  =  16/^  X  sinc(2/,  f)  sinc[2/,(t)-(-T,] 

I~  -  <x. 

Xsinc[2/,(/-|-r2)  ]sinc[2/,(/-)-Tj)  ]  A/ 

and  are  applicable  in  the  same  way  as  the  bispectral  filters. 
The  bi-  or  tri-spectral  masking  filter  is  equivalent  to  a  one¬ 
dimensional  filter  on  only  the  diagonal  factor  in  a  bi-  or 
tri-spectral  approach  to  the  calculation  of  the  bi-  or  tri¬ 
correlation. 


VII.  SAMPLING  REQUIREMENTS  FOR  RECTIFIED 
SIGNALS 

Higher-order  correlation  detectors  and  time  delay  es¬ 
timators  applied  to  certain  signals  have  been  shown  to  ben¬ 
efit  from  rectification  (G.  loup  et  ai,  1990;  Pflug  etai, 
1992b,c).  However,  rectification  can  affect  the  sampling 
requirements  needed  to  avoid  aliasing. 

The  process  of  rectification  means  simply  replacing 
each  negative  amplitude  in  a  given  sequence  with  its  abso¬ 
lute  value,  and  the  unrectified  and  rectified  signal  Fourier 
transforms  have  a  straightforward  relationship.  If  we  sep¬ 
arate  the  discrete  real  sequence  x(t)  into  positive  and  neg¬ 
ative  sequences,  defining  JC+(t)  and  x_(t)  as 

x(t),  if  x:(t)>0, 

0,  if  x(t)  <0 
and 

0,  if  x(t)^0, 

X  ( t)  = 

x{t),  if  x{t)  <0, 

then  the  Fourier  transform  of  x(r)  is  equal  to  the  sum  of 
the  individual  positive  and  negative  transforms, 

and  the  energy  spectrum  is 

ESAf)  =  [X^  if)  +XAf)  ]  [X^  if)  +X  (/)  ]♦ 

=  !^+(/)|'+|^_(/)|'+^+(/)^(/) 


+Xl(f)X_{f). 

If  jc(/)  is  rectified  (denote  the  rectified  sequence  by  x'^) 
then 


x^(t^  = 


x{t),  if  x(r)>0, 

— .x(r),  if  x(/)  <0, 

and  the  transform  of  the  rectified  x(t)  is  simply 
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and  its  energy  spectrum  is 

ES^if)  --  [X^  (/) -X_  if)  1 IX  ^  (/)  -X_  (/)]* 

=  l^+(/)|'+  \XA/)\^-X^(/)X*_(f) 
-X*{f)XAf). 

Since  the  zero-lag  autocorrelation  is  the  same  for  the  un- 
rectihed  and  rectified  signals  and  using  the  central 
ordinate-definite  integral  (sum)  equivalence  (Bracewell, 
1986) 

00  X 

I  £'5,(/)A/=  I  £5f(/)A/ 

y  =  —  X  /=  —  X 

X 

=  I 

/=  -  X 

+  I^_(/)IW 

and 

X 

I  {X^if)X*if)+Xl(f)X_{fmf=0, 

f=  -  CC 

or  equivalently, 

X 

I  {Re[^+(/)]2  +  Re[^_(/)]^  +  Im[^+(/)]2 

/=  -  X 

+  Im[^_(/)]W=0, 

and 

X 

I  {Re[^+(/)]Re[X_(/)] 

/=  -  X 

+  Im  (/)  ]  Im  [^_  (/)  ] }  A/=0. 

Although  this  summation  is  equal  to  zero,  the  individual 
terms  will  not  be,  as  they  are  what  defines  the  difference  in 
the  Fourier  transform  of  the  unrectihed  and  rectified  sig¬ 
nals. 

The  rectification  of  the  whale  transient  shown  in  Fig.  6 
results  in  the  time  series  and  Fourier  magnitude  spectrum 
shown  in  Fig.  12.  Whereas  the  unrectified  whale  transient 
had  /,=32  Hz,  the  rectified  whale  transient  shows  a  top 
frequency  of  approximately  /f  =  150  Hz.  Thus  the  rectified 
autobicorrelation  should  require  about  450  sample  points 
to  avoid  aliasing.  Indeed,  we  see  in  Table  I  that  the  volume 
beneath  the  bispectrum  becomes  constant  to  five  significant 
digits  when  the  signal  is  sampled  with  between  300  and  400 
points.  Similarly,  the  rectified  pulse  contains  higher  fre¬ 
quencies  than  the  unrectified  signal  with  /f=1000  Hz, 
approximately,  and  should  require  about  3000  sample 
points  to  prevent  aliasing  in  the  autobicorrelation.  As 


FIG.  12.  Rectified  whale  transient  and  its  Fourier  magnitude  spectrum. 

shown  in  Table  II  to  five  significant  digits,  the  volume 
beneath  the  bispectrum  becomes  constant  if  a  minimum  of 
3000  sample  points  is  used. 

As  stated,  the  area  beneath  the  spectrum  is  equal  to 
the  central  ordinate  value  of  the  corresponding  correlation, 
and  the  central  ordinate  value  of  the  autocorrelation  of  a 
rectified  signal  is  the  same  as  the  value  when  the  signal  is 
not  rectified.  The  autotricorrelation  of  a  rectified  signal 
also  exhibits  this  property,  consequently  these  values  are 
not  shown  in  Tables  I  and  II  since  they  are  repetitive. 
Therefore,  the  area  and  hypervolume  beneath  the  energy 
spectrum  and  trispectrum  of  rectified  signals  are  not  ap¬ 
propriate  measures  of  the  existence  of  aliasing.  For  in¬ 
stance,  the  rectified  energy  spectrum  of  the  whale  transient 
obviously  contains  aliasing  when  sampled  at  /n=40  Hz 
while  the  unrectified  signal  does  not  (compare  Figs.  6  and 
12  and  see  Table  II).  However,  the  areas  beneath  the  two 
energy  spectra  are  equal.  The  requirements  to  avoid  alias¬ 
ing  in  general  in  the  rectified  autocorrelation  and  autotri¬ 
correlation  are  /N>/f  and  /N>3/f,  ^paralleling  the  re¬ 
quirements  for  the  unrectified  correlations. 

Vill.  SUMMARY 

While  critical  sampling  is  sufficient  to  avoid  aliasing 
when  sampling  the  continuous-time  autocorrelation,  auto¬ 
bicorrelation,  and  autotricorrelation,  avoidance  of  aliasing 
in  calculations  of  higher  order  correlations  and  spectra 
from  discrete-time  data  requires  sampling  rates  higher  than 
the  critical  rate  for  a  bandlimited  energy  signal.  In  these 
cases,  for  the  autobicorrelation,  the  sampling  interval  must 
be  2/3  the  interval  used  for  critical  sampling,  and  for  the 
autotricorrelation,  the  sampling  interval  must  be  1/2  ihe 
critical  sampling  interval.  Examples  have  been  shown  in 
tables  for  both  a  broadband  and  a  narrow-band  bandlim¬ 
ited  transient  using  the  area  or  volume  beneath  the  energy 
spectrum,  bispectrum,  and  trispectrum  (the  correlation 
central  ordinate)  as  a  measure  of  the  presence  of  aliasing. 
This  measure  is  also  used  to  exhibit  the  effective- 
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ness  of  a  bispectral  masking  filter  on  two  sample  signals. 
Using  rectification  as  part  of  a  detection  or  time  delay 
estimation  scheme  for  discrete  data  requires  that  the  sam¬ 
pling  rate  be  selected  or  interpolated  according  to  the  rec¬ 
tified  bandlimit  before  calculating  ordinary  or  higher-order 
correlations. 
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APPENDIX  A:  BISPECTRUM  FOR  A  BICORRELATION 
CALCULATED  FROM  A  DISCRETE-TIME  ENERGY 
SIGNAL 

As  discussed  in  the  body  of  this  paper,  the  discrete¬ 
time  autobicorrelation  could  arise  from  sampling  the 
continuous-time  autobicorrelation  calculated  from  a  real 
continuous-time  energy  signal 


^Ct(1‘1.’‘2)=  JCcT(0-*CT(f  +  Tl)-«CT('  +  T2)‘3?'- 


The  properties  of  the  corresponding  bispectrum  are  dis¬ 
cussed  by  Nielsen  ( 1992)  and  in  this  paper. 

An  alternate  procedure  which  produces  a  discrete-time 
autobicorrelation  is  its  calculation  directly  from  a  discrete¬ 
time  energy  signal.  For  this  purpose  we  make  the  straight¬ 
forward  substitution  of  a  summation  for  the  integral  in  the 
previous  expression 


i>(T|,T2)=  S  x(/)x(r-|-T,)x(/-|-T2)A/, 


with  t,  T],  and  T2  discrete-time  variables  sampled  at  A/. 
Time  domain  functions  not  labeled  CT  in  this  Appendix 
are  discrete-time  functions. 

It  is  important  to  remember  that  the  Fourier  trans¬ 
form  of  the  discrete-time  signal  is  given  by 

with  2tcT(/)  DX(,f),  k/^t=2kff^,  and  X{f)  the  discrete¬ 
time  Fourier  transform  (DTFT)  (Marple,  1987).  The 
DTFT  for  6(t,  ,72)  is 


fc  <c  !  ec  \ 

^(/i./2)=  2  2  f  X  Jc(/)Jc(t  +  r,)x:(/-l-T2)Anexp[-/2iT(/,T,-|-/2T2)]AT,  Ar 

r|  =  -X  72=-*  \  f=-cc  / 


«  /  cc 

X  x:(0:x(f-l-T2)A/ Ar2exp(— /2jr/2T2)  j  X 

71=  —  X  \ 


X  ( t  -I-  r , )  exp  ( -  /27r/,  T, )  At, 


X  X 

=  X  X  x{t)x(.t-\-T2)i^t^T2txp{—i2nf2T2)X{f^)e\p(+i2vf\t) 

f=-*x  r2=— X 


X 

=  X  Jc(t)Arexp(-i-/2ir/|r) 

f—  —  X 


X  x(f-|-T2)exp(-/27r/2T2)AT2  p( /,) 

T2=-*  / 


X 

=  X  Jr(OA/exp(-t-/2jr/,f)[A^(/2)exp(-(-i2ir/2/)]^( /i) 

/=  -  X 


X 

=  X  xU)tXp[+aiT(f^+f2)t]MX{f0X(f2) 

•  =  —  X 

and  for  x{i)  real 


X  x(/)exp[ -/27r( /,-|-/2)t]Ar  X{  f^)X{  f2) 
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^X*(f^+f^)X{f^)X{h) 

X(  fi  +  Iki  /n  )  j 

xf  2  /i+/2  +  2^3 /n)  y 

\*3=-0=  J 

It  is  the  transform  replications  in  this  expression  that  pro¬ 
duce  the  aliasing  in  the  TAD  of  the  transform  of  the 
discrete-time  autobicorrelation  calculated  from  discrete¬ 
time  energy  signals  as  described  in  this  paper. 


APPENDIX  B:  BISPECTRUM  FOR  A  BICORRELATION 
CALCULATED  FROM  A  DISCRETE  TIME 
STOCHASTIC  POWER  SIGNAL 


A  result  corresponding  to  that  of  Appendix  A  can  be 
derived  for  stationary  stochastic  power  signals  following 
the  derivation  of  Hinich  (1990)  [his  Eq.  (3.4)  and  subse¬ 
quent  discussion].  Hinich’s  derivation  is  for  continuous- 
and  discrete-time  autobicorrelations  (equivalent  to  autoco¬ 
variances  for  his  zero  mean  assumption)  calculated  from 
continuous-time  stationary  stochastic  signals  and  for  their 
spectra.  Hinich’s  Eq.  (3.4)  is 

E{AA  f)A,(g)A,m 


=E  JJJ x(/)x(r)jc(i)exp[-/'2tr(  ft+gr 


+  hs)  ]dt  drds 


=  JJ C;f;j^(r,0)exp[-/27r(  fT+gd)]dT dd 

X  {  en-pl-ilvi  f+g+h)t]dt 


=  B,{f,g)8(f+g+h), 
for  which 


A,c=  I"  x{t)e\p(—i2nft)dt 

(A^  a  complex  random  variable), 

=E{x(t)x(t  +  T)x(f+0)}, 

and  B^{f,g)  is  the  autobispectrum  of  x{t).  Hinich  goes 
on  to  say,  “It  follows  from  (3.4)  that 
/  B^(f,g)df  dg=E{fA^(f)A^(g)A^(,h)dfdg}  when 

h=—f—g.  But  since  /4^(/)=0  for  />/o,  then 
*  Bj,if,g)=0  for  /-|-g> /o”  He  then  discusses  the  conse¬ 
quences  of  sampling  the  continuous-time  The 

^  results  correspond  in  terms  of  aliasing  to  those  given  in  this 
paper  for  sampling  the  continuous-time  autobicorrelation 
of  energy  signals,  and  agree  with  those  of  Nielson  ( 1992) 
as  well. 


If  instead  an  analogous  result  is  desired  for  discrete¬ 
time  autobicorrelations  calculated  from  discrete-time  sta¬ 
tionary  stochastic  power  signals,  the  corresponding  deriva¬ 
tion  is  as  follows.  Let  Xqx(^)  the  discrete-time  stochastic 
signal,  and  use  a  tilde  to  represent  quantities  calculated 
from  it.  [Contrary  to  Appendix  A,  we  let  the  unsubscripted 
x(/)  be  the  continuous-time  function  to  agree  with  Hin¬ 
ich’s  notation.]  Then 


X 

XDT(t)DAA/)=  1  A,{f-lkf^), 

k=~cc 


and 


?«;.(  ^  )-*Dt(  f  f  +  ^)  }■ 

Therefore 

E{A,{f)AMUAh)} 

CX  X  X 

—  X  —  X  T=  —  X 

Xexp[— i2ir(  ft-^gr-\-hs)\^t  Ar  Ar] 


X  X 


=  X  X  c^(T,d)txp[-n.Tr{  fr+gd)] 

Ts=  —  X  0=  —  X 
X 

XAtA0  X  exp[  — /2ir( /-l-g-f-/i)r]Ar 

f=  —  X 

=  B,(/,g)6(/-fg4-A). 

From  this  result,  B^(f,g)  corresponds  to  E{A^{/ ) 
XA^ig)A,(-/-g)}.  But 


AAf)^  X  AAf-2k,f^), 


1--X 


AAg)=  X  AAg—^kiffi), 


and 


^xi-f-g)=  X  AA-f-g-2kiffi). 

^3=  —  X 

One  can  no  longer  assert  that  since  A  A  —  f —g)  is  zero  in 
the  domain  outside  (/+g)  >  /n,  the  autobispectrum  will 
be  zero  in  that  domain  for  the  autobicorrelation  calculated 
from  discrete-time  stationary  stochastic  power  signals.  In¬ 
stead  one  must  look  at  A^(/—g),  which  is  not  zero  for 
(/-rg)  >  /n.  Therefore  the  PAP  (the  OT  of  Hinich)  and 
TAD  will  be  nonzero  for  this  bicorrelation  even  for  sta¬ 
tionary  signals,  if  A/=l/(2/,),  with  /,  the  highest  fre¬ 
quency  present  in  x(r).  Only  if  Ar>l/(3/,)  will  the  PAP 
and  TAD  be  zero.  This  result  is  in  contrast  to  that  for 
sampling  the  continuous-time  autobicorrelation,  and  the 
corresponding  Hinich-Wolinsky  (1988)  test  for  stationar- 
ity,  for  which  the  PAP  (OT)  is  zero  if  the  autobicorrela¬ 
tion  is  sampled  at  A/=  1/(2/,). 
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